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Abstract 

The renormalization group is used to resum leading logarithmic contributions 
of the form a" +1 /?o log n (A///) to the gap equation appropriate for high den- 
sity QCD. The scale dependence of the strong coupling constant a s increases 
the gap by a factor of exp — 1^ ~ 20 in the small coupling limit. 
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The present efforts towards understanding QCD at extremely high density |T[- |T(J are 
an important step toward understanding QCD at moderate densities which may one day be 
accessible to experimental study. For particular combinations of quark colors and flavors 
at high densities, it is likely that a superconducting gap breaks color and flavor symmetries 
in interesting ways. Although this symmetry breaking is nonperturbative, it occurs when 
QCD is weakly coupled, and therefore perturbative QCD (pQCD) can be used to derive 
properties of the superconducting phase. The most basic property is the size of the gap 
itself. In beautiful work, Son || realized that the gap in pQCD is dominated by magnetic 
gluon exchange, regulated in the infrared by dynamic screening. He arrived at an estimate of 
the gap by renormalization group (RG) arguments, and found A ~ cg~ 5 exp -j^j-j in the 
small coupling limit, where g s is the strong coupling constant evaluated at a scale of order 
the chemical potential, fi, and c is a numerical coefficient. This surprising result was also 
obtained directly from the Schwinger-Dyson (SD) equations ||. Further, numerical solution 
of the gap equation || with particular boundary conditions yields values of the gap that are 
within an order of magnitude of the small coupling result for the same boundary conditions. 

In order to refine the perturbative estimate of the gap, the coefficient c must be deter- 
mined. At this one-loop order in the SD equations there are several important contributions 



to c [e.g. |Tl| ], one being the ambiguity in the scale at which the strong coupling constant 
g s is evaluated. Working in an effective theory appropriate for momentum transfers below 
2/i, it has been suggested that choosing a renormalization scale A = \x (mid way between 2fi 
and 0) provides an estimate of the contribution to the running of g s ||. While plausible, 
in order to determine c a somewhat more rigorous estimate is required, and in fact we will 
show that changing the scale of g s in the solution of the gap equation does not capture the 
largest contribution. The scale ambiguity in the gap equation, as discussed in will 
be reduced if computations are performed at higher orders in the g s expansion, i.e. 0(a 2 s ). 
For many processes, an estimate of the scale can be made by computing terms of the form 
a 2 /3o, {(3q = 11 — |n/ in vacuum, rif being the number of dynamical quark flavors) from 
the gluon vacuum polarization diagrams arising at two-loops [p~|]. Introduced by Brodsky, 
Lepage and Mackenzie (BLM), this scale setting procedure has been used extensively for 
several different processes, including the inclusive decay of heavy quarks, e.g. [13|]. However, 
there is a problem applying this technique to the gap equation because log A is an expan- 
sion in g s and not a s . In this work we will compute and resum contributions of the form 
a n+1 p^ log n (A//i) to the gap equation for A using the renormalization group. 

Explicit construction of the effective field theory with which we work, appropriate for 
momentum scales below 2/i, can be found in papers by Hong ||. Neglecting higher order 
coupling between the SD equations relating the gluon and quark two-point functions, only 
the graphs shown in Fig. (|l|) need be calculated in order to resum the leading logarith- 
mic contributions. Regularizing the perpendicular gluon momentum kj_ with dimensional 
regularization and minimal subtraction one finds that at A ~ 2/x the gap equation is 
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We have neglected the (po ~ Qof terms in the gluon propagators as they constitute higher 
order corrections in the effective field theory. In order to obtain eq. ([[]), we have used the 
small go limit of the Landau damping term, 
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The "constant" C is the coefficient of a four-quark operator in the effective theory whose value 
is determined by matching QCD to the effective field theory. In order to absorb the scale 
dependence of the one-loop graphs, C has precisely equal and opposite scale dependence. 
At higher orders in the and g s expansion four-quark operators will enter with scale 
and gauge dependence that precisely compensates the scale and gauge dependence of the 
loop graphs in the effective theory. With no further mention of the gauge dependence of 
C, we exploit the scale dependence of C to resum all leading log contributions to the gap 
equation. In what follows we neglect the scale dependence of the g s 's that appear in eq. (0), 
and assume that they are evaluated at 2/i, as variations from these terms are higher order 
effects. In fact, it is simpler to separately evolve the contribution from electric and magnetic 
gluons by setting C = C^ E ' + C^ M \ The assertion that A is scale independent leads to a 
differential equation for the 
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where 7® = 2 for both % = E, M. Straightforward solution of this differential equation leads 
to a RG improved gap equation, where we chose A = A# for the electric and A = Am for the 
magnetic gluons) 
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An important point to note here is that the gap equation does depend upon the chemical 
potential /x, however, it is independent of the renormalization scale A in the effective field 
theory. While A(p ) is symmetric under p — ► — p , the integrand is not symmetric under 
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Qo ~Qo- In the region where both electric and magnetic gluons are effectively massless, 
\e < A < 2/i, the evolution of the strong coupling is determined by /3 = 11, which can 
be seen by considering the gluon three-point function and making use of gauge invariance. 
However, the /3-function in the region below A#, where the electric gluons are not dynamical, 
is yet to be computed; in eq. (^) and throughout this work we will assume that the magnetic 
and eletric /?o's are the same. 

Proceeding along lines similar to Son ||, a second order differential equation for A(po) 
can be obtained, and solved for given boundary conditions. However, in contrast to Son 
who sets up a differential equation in the variable log (po), we use the variable a s (\% 3 po ), 
which at leading order recovers the Son result, and enables us to include the higher order 
contributions from the evolution of g s . 

Assuming that go ^> A, we can write the gap equation in eq. (W) as 
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where au = a s (A^ 3 |fco| 1//3 )- Using the somewhat brutal approximation of ||, 
oi s {\po ~ qo\) = 9(a s (q ) - a s (p ))a s (p ) + 9(a s (p ) - a s (q ))a s (q ) 
we are able to convert the integral equation into the differential equation 
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The general solution to this differential equation is 
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where Jq(x) and Y Q (x) are Bessel and Neumann functions, respectively. In the limit where 
a a is small, a self consistent solution where g s log(A//i) ~ 1 can be found. Taking the 
asymptotic limit of the Bessel functions in eq. (H), 
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where A and (p are constants. As expected this reduces to the solution found by Son when 
the argument of the Sine function is expanded in powers of logpo, but with the advantage 
that we have been able to resum the leading logs directly. 

In order to determine A and (ft we return to the integral equation in eq. (|5|). Given the 
form of the gap in eq. @, it is clear that the integral on the right hand side of eq. (§) is 
divergent, as a — > 0, for any value of 0, and hence needs to be regulated in the ultra-violet. 
While unrelated, similar integral equations appear in effective field theory studies of the 
three-body problem fll4} . Efforts to dimensionally regulate such integrals have failed up to 



4 



this point, however, cut-off regulation has been used successfully |T^] . Therefore, we will 
introduce a cut-off, a a, on the longitudinal direction to define the integral equation and 
consequently, the countertern C(2/z) becomes dependent on this cutoff in the longitudinal 
direction C(2/i) = C(2fx, a a) Using the approximation in eq. (@) and 7^ = 7^ = 2, eq. @ 
becomes (including the infrared cut-off defined by A, a^) 
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It is clear from eq. (|l0|) , that A(a) = 0, when a a = «■ As the result is required to be 
independent of a a, we find that, in the small a a limit 
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It is important to notice the dependence of C(2/x, qa) on the oa- Unlike the counterterms 
in perturbative effective field theories, this counterterm changes rapidly with scale (except 
near oia = a?), and is periodic, as found in three-body physics [Q. In the small coupling 
limit, C(2(i,a A ) ~ C(2/x,a) ~ # 2 (2/x). 

As the solution to the gap equation is independent of oa, by construction, we choose to 
work at a a = a. We find 
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From eq. ( |T0D it is clear that a p = a a when jjA. = o ; which leads to 
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is expected to be a number of order unity. The first two factors constitute the result obtained 



by § and §. The third factor exp (^J 
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20 arises from the scale dependence 



of the strong coupling constant in addition to the logarithmic dependence of the one-loop 
gluon graphs themselves. Somewhat surprisingly, it is independent of g s . This factor will 
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not be recovered by simply evaluating the leading term at a renormalization point between 
and 2/i, as such corrections are of the form exp (rg s log (A//z)), where r is some number. 
The fourth factor is of order unity, and must be determined by matching to a full QCD 
calculation, which we have not done. 

To make sure the enhancement found above provides a good estimate of the true effect of 
the running of g s we have numerically solved eq. (||) for a a ~ «• As the full QCD calculation 
has not been performed we set C(2/x,2/z) ~ C(2/i, a^) = 0, since it is not enhanced by j3 . 
For (j, = 10 6 MeV (n f = 2) and a A = a we find A = 4.4 MeV for p = 0, and A = 32.5 MeV 
for (3q = 11. For this density A is enhanced by a factor of ~ 7 by the running of g s . Similarly, 
for fi = 10 4 MeV and a A = a we find A = 5.7 MeV for (3 = 0, and A = 43.0 MeV for 
(3o = 11. The difference between the estimated enhancement (~ 20) and the one found 
numerically (~ 7) may be attributed to the approximation in eq.(^|). It is reassuring that 
A changes by only ~ 2 when the cut-off momentum is increased by an order of magnitude, 
while keeping C(2/i, 2/x) ~ C(2[i, a a) = 0. 

In producing these numbers the go ^ |q| limit of the gluon screening functions was 
employed. While we expect that such an approximation will modify the results by a factor 
of order unity, we do not expect the enhancement due to the running of g s to be substantially 
different. 

In conclusion, we have found a significant enhancement in the size of the color- 
superconducting gap, A, due to the scale dependence of the strong coupling constant, g s . 
We have resumed leading logarithms using the renormalization group and obtained a scale 
independent gap equation to leading order in the effective field theory expansion in g s and 
l//i. Terms that do not involve the large factor of (3q have not been calculated in this work, 
but are naively estimated to be smaller. 

Matching to the effective theory needs to be performed at higher orders to make a precise 
prediction for A. In our calculations, we have assumed that both the magnetic and electric 
gluons have the same /3-function for all scales below 2/z. Below the Debye mass, the (3- 
function for the magnetic gluons will deviate due to the absence of logarithmic terms from 
electric gluon loops. However, the difference is expected to be of order 20% and our work 
should be taken to be only an estimate of this effect. 
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